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Abstract 

We give a comprehensive analysis of how scalar and tensor perturbations evolve in cosmologies 
with a smooth transition from power-law-like and de Sitter-like inflation to a radiation era. Analytic 
forms for the super-horizon and sub-horizon perturbations in the inflationary and radiation dominated 
eras are found. 

I Introduction 

In Jl| , the evolution of perturbations for a decaying vacuum cosmology given in Q , with smooth exit 
from inflation to radiation-domination, was found exactly on super-horizon scales. The solution confirmed 
the standard result §, |, §, g, § within gauge invariant perturbation theory that super-horizon density 
perturbations are strongly amplified through the transition from inflation to radiation. This transition is 
often approximated as an instantaneous jump Q 0, since super- horizon modes change on a timescale 
that is much greater than the transition time. To avoid problems in the matching conditions for a jump 
transition, we consider only smooth transitions. 

Given the variety of methods and models that have been used, we believe it is useful to do a com- 
parative study, using the same method on two broad classes of transition models, i.e. those based on 
power-law-like and de Sitter-like inflation. We use the class of power-law-like models first given in 
The second class of models is a generalisation of the models of || |) with a de Sitter-like inflationary 
era. The models, together with the neccessary equations in perturbation theory, are discussed in sec- 
tion 2. In section 3 we present our results for scalar and tensor perturbations on super- horizon scales, 
and show how they confirm the standard predictions. In section 4 we present solutions for scalar and 
tensor perturbations on sub-horizon scales. However, toy models that evolve smoothly from inflation 
to radiation dominated eras are only applicable for super-horizon modes around the time of transition, 
since on sub-horizon scales, the dynamics of the reheating era have a significant effect. The solutions 
for sub-horizon scales therefore do not apply during the transition era. Section 5 gives some concluding 
remarks. 

Throughout this paper we use units in which 8irG = 1 = c. 



II The models and perturbation equations 

A flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background has metric 

ds 2 = -dt 2 + a 2 {t)[dx 2 + dy 2 + dz 2 ] , (1) 

where t is the physical time and a(t) is the expansion scale factor. The background pressure p, energy 
density p and Hubble parameter H = a/ a (where a dot denotes a derivative with respect to cosmic time 
t) satisfy the field equations 

p = -2H-3H 2 , (2) 
p = 3H 2 . (3) 
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The pressure index and effective adiabatic sound speed are defined by 
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Inflation is characterised by d > 0, so that the epoch t c of exit is given by d(t c ) = 0. The radiation 
dominated era is characterised by p « lp, H rs l/2t. 

Scalar and tensor perturbations lead to a metric of the form 



ds 1 



a 2 (?7){-(l + 2$)ffr/ 2 + [(1 - 2*) 7y - + hi^dx 3 } . 



(6) 



in the longitudinal gauge, where r\ is conformal time (dr\ = dt/a) Q. Here $ and \I/ are gauge- invariant 
amplitudes of the scalar metric perturbations. The field equations require <£> = ^ when the spatial part 
of the energy-momentum tensor is diagonal. The tensor perturbations are gauge-invariantly described 
by hij — hQij, where Qij is a symmetric traceless three tensor and h represents the amplitude of a 
gravitational wave. 

The evolution of adiabatic scalar perturbations $ is given by B 



rf 2 $ 
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3aH{l + c 2 ) 



d<P 

drj 



\7 2 <f> + 3a 2 H 2 (c 2 s - w)<5> = 0. 
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Decomposing $(77, x) into eigenmodes $(77, k) of the comoving Laplacian, and then using the cosmic time 
t and scale factor a as the dynamical variables respectively, gives 
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where fc is the comoving wave number. The decomposition allows us to follow the evolution of a single 
mode through the inflationary and radiation dominated eras. The physical length scale A correponding 
to the comoving wavenumber k is A = 2ira/k. This will equal the Hubble radius H~ l when 



k = 2iraH . 



(10) 



Super-horizon modes are characterised by k <C aH. For these modes one can neglect the k 2 term in 
equations (||) and (||) (and subsequent equations). The solutions of the truncated equations will be a 
good approximation to the limit of the solutions of the full equations, provided pt] 



fc 2 $ -> Oas k -> 0. 



(11) 



For adiabatic super-horizon scalar perturbations, the growing modes have a conserved quantity pp| which 
can be used to determine $: 




H f aH 

Hp 



dt. 
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Equation ( |l2|) can be used to express the growing perturbations at late times in terms of their early 
time forms [Bj . In || [| , £ is used to estimate the amplification of perturbations through a smooth 
transition. This amplification has been disputed in [jllj (see also |L2|, H). 
The density perturbations in the longitudinal gauge are given by pf: 



The Fourier mode h of the amplitude h of the tensor perturbations satisfies the equation 

d 2 h „ „dh ( fc x " 



which may also be re-written with a as the dynamical variable: 



2 d 2 h ( a 2 dH\ dh k 2 - 
a ^ + { Aa+ H^)d-a + ^p h = °- (15) 

For super-horizon modes one can reduce to quadrature: 

I-o+dJS-o+dJ*,. (.6, 

where C and D are constant, ft follows that D corresponds to the decaying mode and therefore there is 
no amplification in long- wavelength tensor perturbations |ij . 

A Power-law-like models 

In (s| a simple class of power-law-like inflation models, which exit smoothly to a radiation era, is given 

by 

a (t) = -£=(1 + m)( 2 »- 1 )/ 2 t"(i + mte)^- 2 ")/ 2 , (17) 
vis 



where the parameter n satisfies n > 1 to acheive inflationary expansion, and 



n + Jn(2n - 1) 

771 = 

2n(n-l) 

For t <C mt e , 

/ 1 \ (2n-l)/2 / ^ 

showing the early power-law-like inflation era. For t ^> mi c , 

as in a radiation dominated era. Thus we have a cosmology which evolves smoothly from power-law-like 
inflation to a radiation dominated era. The Hubble rate follows on differentiating equation (|l7j): 

™ = M' (18) 



Using equations (|^) and (||), equation (|T^) gives 

_ t 2 + Anmtet - 4n(3n - 2)m 2 t 2 

P ~ 4t 2 (t + mi ) 2 ' ( } 

3(t + 2nmt c ) 2 , , 

p = WW' (20) 

It follows that the pressure index and effective adiabatic sound speed given by equations (Q) and (|5|) are 

t 2 + Anmt c t - 4n(3n - 2)m 2 t 2 e ,„„ . 

w = ~ o — o — s— , 21) 

3i 2 + 12nmt e t + I2n 2 m 2 t c 2 

2 _ i(£ + mt c )(t + 2nmt c ) — [t 2 + 4nmt c t — 4n(3n — 2)m 2 t 2 ](2t + mt c ) 

Cs ~ 3t(t + mt c )(t + 2nmt c ) ~ 3(i + 2nmt c ) 2 (2t + mt c ) ' ^ ' 
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The evolution of scalar and tensor perturbations for this model is given by equations (Q), ( |13[) and 
(pi). Using equations (§]]), @, @ and (||), we get 



where 



and 



and 
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rf<E> 



t + 2nmt c 
2{t + mt e ) 



V(i)$ = 0, 



5t(t + mi e )(i + 2nmt e ) - (5t 2 + 20nmt c t + 4n(n + 2)m 2 t 2 c )(2t + mt c ) 

t(t + mt e )(t + 2nmt c ) - (t + 2nmt c ) 2 (2t + mt e ) 

t(t + mt c )(t + 2nmt c ) - [t 2 + Anmt c t - An('in - 2)m 2 t 2 ](2t + mt e ) 

t(t + mt c )(t + 2nmt c ) - (t + 2nmt c ) 2 (2t + mt c ) 

t 2 + 4nmt c t - 4n(3n - 2)m 2 t 2 c 2 i c (l + m) 1 - 2 ™* 2 " 2 "^ + mt c ) 2n+1 
h K 



— = -2 



(t + 2nmt e ) 2 

4fc 2 t e (l + mf- 2n t 2 - 2n {t + mt c ) 2 " +1 l 
3a 2 (t + 2nmt c ) 2 + 



a 2 (t + 2nmt c ) 2 



4t(t + mt e ) 



t + 2nmt c 



~dt 



(23) 



(24) 
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2n— 1+2— 2ra1 



ft = 0. 



(25) 



B de Sitter-like models 

In H a new exact inflationary solution is given, using a as the effective dynamical variable. The Hubble 
rate is given by 

H(a) = ^r4, (26) 

where £f c is the Hubble rate at exit. A similar model [but one which does not produce an exact solution 
for super-horizon scalar perturbations, unlike ( |2^ ) pi) is proposed in S: 

ol/2 rr „2 

ff (°) - r 4 1 4°W2 ■ ( 2? ) 

(a 4 + a 4 ) 1 /^ 

For both Hubble rates, when a <C a c , H is approximately constant as in a de Sitter inflation era, and for 
a a c , the Hubble rate decays like a -2 , as in a radiation dominated era. Thus for both models we have 
a cosmology which evolves smoothly from inflation to a radiation dominated era. 
Here we generalise these models by considering the class of Hubble rates 

2 1 l t a 2 H 



parametrised by I (I > 0). For equation (M), I = 1, and for (J27|), 1 — 2. Using equation (|28|), H ~ 
constant for a <C a c (inflation), and -ff oc a^ 2 for a 3> a (radiation domination), confirming that the 
model has the desired properties for all positive I. 
Using equations (||) and (||), equation ( pS| ) gives 

o2/^ fr2„4/'„2£ o_ 2£\ 

P (a^ + a 2 1+2 / f ' 1 ' 

2^H 2 aj 

9 (a 2 ? + a 21 ) 2 ' 1 ' [ ' 
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The pressure index and effective adiabatic sound speed given by equations (0) and (|^) are 

1 fa 2t -'iaT 



a M - (21 + 3)o, 



2t 



(31) 
(32) 



Applying equation di"o| ) to equation fl2S| ) and taking fc c = 2ira c H c (the comoving wavenumber of the 
Hubble radius at exit), the epoch a_ of leaving and a + of re-entering the Hubble radius are given exactly 
in this class of models by 



a± _ k c 
a e k 



2( 



(33) 



Scales with k > k c never cross the Hubble radius and remain sub-horizon. For super-horizon scales with 
k <C fe e j equation (|33|) reduces to 



a± 



2V^ f 



±i 



The scalar and tensor perturbations for this model are governed by equations 
Using equations (28), (|3l| ) and ( j32| ) we get 



© and ©. 



4a ; 



2C 



(2£-2)a, 



2/ 1 
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k\ 2 (a 2e + a 2 J) 2 / e 2£a. 
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2 2 / e a 2 a 2 
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(a M + ) 

2 2/£-2 fl 2 a 2 



h = 0. 



(34) 



(35) 



(36) 



III Super-horizon perturbations 

For the modes which leave the Hubble radius (during inflation), while they remain outside the Hubble 



radius, we can neglect the /c-term in equations ([23|), (|2J|), (|25|), (|34|), ( |35D and (|36|). Only equations fl25|) 
and (|36|), using ([if]), give exact solutions for the truncated equations. However, using jl4. |l5) we can 
give analytic forms for the solutions in the inflationary and radiation dominated eras for equations (p3|), 
(pif), ( |34| ) and (j35|), since they reduce to Bessel forms. For all solutions, Aj, Bj, j = 1,2, ... are used to 
denote arbitrary constants. For equation ( [Tl| ) to be satisfied it is required that for scalar perturbations 
Aj, Bj ~ k q as k — ► 0, where q > —2. For all scalar perturbation solutions this is assumed to be the case. 

A Power-law-like model 

For t <C mt c equation ( |23| ) has solution 

, \ n+l 



$ ~ Al + Bl [ -| ) , (37) 



describing modes that leave well before exit, while for t 3> mt e 

N 3/2 



$ » A 2 + B 2 ( -f ) , (38) 
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describing modes that remain super- horizon in the radiation era. In equations ( |37| ) and (|38|), Aj corre- 
spond to the constant modes, and Bj to the decaying modes. 

The conserved quantity given by equation (|l2|), applied to the case of the model given by (|l7|), leads 

to 

~_ (t + mt c )( 2n -V/ 2 (t + 2nmt c ) f t n (t 2 + 4nmt c t + 2nm 2 t c 2 ) 
~ C ¥+ l J (t + 2nmt e ) 2 (t + mi )( 2 ™- 1 )/ 2 ' 

where we neglect the decaying contributions (and therefore no constant of integration arises from the 
integral) . Amplification of the scalar growing modes according to equations ( |37| ) and ( |38| ) is given by the 
factor a = A2 jA\ . Using equation (|3^) , we find that 

M = C 



(n + 1) ' 
and 

Therefore 

a=|(n+l)>|, (40) 

so that amplification is only significant for large n. This is in agreement with results based on instanta- 
neous transition models (|, |[ Q . 

The density perturbations in the inflationary and radiation dominated eras, are found analytically, 
using equations (^4|), (|37]) and ([38]). For modes that leave the Hubble radius well before exit (t <C mt e ) 

^-2A 1 + WM" + \ (41) 
P n V t ) 

while for modes that are still super-horizon during radiation domination (t 3> mt e ) 

S -P^-2A 2 + 4B 2 ( t -2-) 3/2 . (42) 



P \t 

The density perturbations undergo the same amplification (^) as the potential perturbations. 
The tensor perturbations are given by equation (|l6|): 

h = A 3 + B3J t- 3n (t + mt e ) 3(2 "- 1)/2 * . (43) 

During power-law inflation (t <C mt e ) 

h*A 3 + B 3 m - * e (±) , (44) 



1 - 3n V t , 

while for those modes that remain super- horizon during radiation domination (t 3> mt c ) 

b A 3 - 2B 3 t~ 1/2 . (45) 
The results (ff3|)-(f45|) imply that tensor perturbations do not grow on super-horizon scales. 

B de Sitter-like model 

For modes that leave the Hubble radius when a <C a c , equation (|34|) has super-horizon solutions 

/ \ 2t 

iw^f-i) + B 4 ^. (46) 
V a c / a 
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For super-horizon modes during radiation domination (a 3> a, 

3/2 



4> 



(7)' 



3/2£ 



Vf(?)' 



(47) 



where Z v denotes a linear combination of the Bessel functions J v and Y v . Using the asymptotic forms of 
the Bessel functions Jul fL5|, equation (E^) reduces to 



(48) 

It follows that in the inflationary era, |$| grows as a 2e , while in the radiation dominated era, |$| is 
approximately constant (while the scales are still super- horizon). Thus growing super- horizon scalar per- 
turbations are strongly amplified during inflation and then remain approximately constant after inflation. 
The amplification factor is clearly much greater than that for the power-law-like model. 
We can calculate the amplification using equation (|lj), which gives 



2C 



a{a 21 + af) 1 ' 



(a 2 * + a 2 /) 



da . 



(49) 



Evaluating the integral in the inflationary and radiation era's and using equations ( ff6[ ) and ( (48[ ) , we find 

2( 



A 4 



and 



-4, 



2£ + 1 
2C 



The amplification from the epoch ai (a; <C a c ) to the radiation era is thus 



1 



a=-(2£+l) ^ 



a; 



(50) 



The explicit dependence in (|50|) of the amplification on the initial epoch aj shows that amplification, and 
thus the spectrum, is not scale invariant. 

The density perturbations are given using equations (|35|), (Q) and ([17]). For modes beyond the Hubble 
radius during inflation 



21 



P \a c 
while for super- horizon modes during the radiation dominated era (a> a c ) 



(51) 
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-*(2l-3)/2l 



2 ta 



- 2 



(7) 



3/2 



-3/21 



which leads to 



Sp 



p v a J 



(52) 



(53) 



During inflation (a <C a c ), <5p//0 grows as a 2f , while in the radiation dominated era (a 3> a c ), <5p//J is 
approximately constant. The amplification is the same as (|50|). 

The evolution of the tensor perturbations for super-horizon modes is given by equation (fHf): 



h = A 6 



B 6 



(a 2e + a 2 c e )^a c 



da . 



(54) 



7 



For the inflationary era (a <C a c ) equation ( p4[) gives 

fc«A 6 -±B 6 (^) 3 , (55) 
3 V a / 

while for the radiation dominated era (a a c ) 

h*A 6 -B e (^) . (56) 
Unlike the scalar the tensor perturbations do not grow during inflation on super-horizon scales. 



IV Sub-horizon perturbations 

For perturbations on scales inside the Hubble radius (where k/(aH) is not negligible), equations (p3|), 
(p4|), (p5|), (p^), ( |35| ) and ([36j) can not be solved exactly. We can however give analytic forms for the 
solutions in inflation and radiation domination using |f4|, |f5|, since the equations reduce to Bessel forms. 

A Power-law-like model 

Defining the comoving wave numbers 



fa = 

k 2 = 



a c (l - n) 



< c (l + l/m)( 1 - 2 ")/ 2 ' 



2t c (l + m)( 1 -2»)/2 ' 

we find the following. During power- law- like inflation (t <C mt e ), the scalar and tensor perturbations 



evolve as 
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(57) 



2(1 -n) fc /* y 1 - 3 ™)/ 2 
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(3n-l)/2(l-n) 



- 2 
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Z, 
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fc / t 



ki V t, 
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ki \t e 



fa \t c 



(58) 
(59) 



During radiation domination (t 3> mt e ), we find the following: 



4> 



3/4 
- I Z: 



3/2 



3/2 



k ft 



k 2 \t c 
Aj sin 
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(61) 
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(62) 



As expected, density perturbations oscillate after re-entry, since radiation pressure balances gravitational 
infall. Note that while there are non-decaying oscillations in the density perturbations, the tensor per- 
turbation oscillations are purely decaying. At late times, the density perturbation (^TJ) has the form 



5_p 
P 



A cos 



fc a 
k 2 a c 



B sin 



fc a 
fc 2 a c 



in agreement with the standard results for the radiation era |ll| (p. 152). 
B de Sitter-like model 

Sub- horizon perturbations during inflation (a <C a D ) have the form 
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Sub- horizon perturbations during radiation domination (a » a c ) are given by 
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(67) 



(68) 



Once again, we see that the density perturbations contain non-decaying oscillations while the tensor 
perturbation oscillations are purely decaying. At late times, the density perturbation (|67j) has the form 



— — ~ A cos 



2(1-*)/* fc c a c 



5 sin 



2(1-*)/* fc c a c 



as in the power-law-like case. 



V Conclusions 

Our main result is a comprehensive presentation of analytic solutions for the scalar and tensor pertur- 
bations on super- and sub-horizon scales for two classes of cosmology, where the scale factor evolves 
smoothly from an inflationary era to a radiation dominated era. These classes encompass parametrised 
power-law-like and de Sitter-like inflationary behaviour. Our results confirm the amplification of super- 
horizon scalar perturbations for all values of the parameters which define the classes of models. Explicit 
forms, given in equations (M) and (pOh , of the amplification for scalar perturbations show how the de 
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Sitter-like inflation is much more effective in amplification than the power-law-like inflation. On sub- 
horizon scales, the detailed form of oscillating modes after re-entry was found, showing that there are 
non-decaying modes for density perturbations, but not for tensor perturbations. Taken together, these 
results provide a unified catalogue of perturbation solutions in what can be taken as the two main types 
of inflationary models with smooth exit. 
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